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Abstract 

Background: The Alberta oil sands are an important economic resource in Canada, but there is growing concern over 
the environmental and health effects as a result of contaminant releases and exposures. Recent studies have shown a 
temporal and spatial trend of increased polycyclic aromatic hydrocarbon (PAH) concentrations in sediments and 
snowpack near the Athabasca oil sands operations (i.e., open pit mines), but thus far similar studies have not been 
done for the Cold Lake region where steam assisted gravity drainage (in situ) extraction is performed. 

Methods: Many PAHs are known mutagenic carcinogens, and this study measured soil and atmospheric 
concentrations of PAHs in the Cold Lake region to assess the excess lifetime cancer risk posed to the First Nations' 
inhabitants of the region. Using both deterministic and probabilistic risk assessment methods, excess lifetime cancer 
risks were calculated for exposures from inhalation or inadvertent soil ingestion. 

Results: The mean excess cancer risk for First Nations' people through ingestion who engage in traditional wilderness 
activities in the Cold Lake region was 0.02 new cases per 100,000 with an upper 95% risk level of 0.07 cases per 100,000. 
Exposure to PAHs via inhalation revealed a maximum excess lifetime cancer risk of less than 0.1 cases per 100,000. 

Conclusions: Excess lifetime risk values below 1 case per 100,000 is generally considered negligible, thus our analyses did 
not demonstrate any significant increases in cancer risks associated with PAH exposures for First Nations people 
inhabiting the Cold Lake region. 
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Background 

Canada has the third largest proven oil reserves in the 
world [1]. The oil sands in Alberta represent the vast ma- 
jority of Canada's oil reserves, and are the largest source of 
crude bitumen in the world. Currently, there is much 
speculation in the media, and within the general public, 
concerning the impact of the Alberta oil sands operations 
on the environment and human health. Some areas of the 
Alberta oil sands have very high background concentra- 
tions of contaminants such as polycyclic aromatic hydro- 
carbons (PAHs) [2-4], but there is increasing evidence that 
there have been significant increases in environmental 
levels of PAHs and metals as a result of the industrial de- 
velopment of the oil sands [5-11]. Additionally, there is 
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concern in the general public and media about increased 
incidence of cancer and other health effects in areas im- 
pacted by the Athabasca oil sands [12,13]. 

There are three major oil sands regions in Alberta: (1) 
Athabasca (the largest); (2) Cold Lake; and (3) Peace 
River. The Cold Lake region, the focus of this study, is 
the largest in-situ thermal heavy oil (bitumen) extraction 
operation in the world, producing 54,600 barrels of bitu- 
men each day. Currently, in-situ bitumen extraction ac- 
counts for 47% of oil sands extraction in Alberta, and 
will soon become the dominant form of oil sands extrac- 
tion, overtaking open pit mining [14,15]. There is 272.3 
billion barrels of known oil reserves in the Alberta oil 
sands, and 93% of the known reserves can be recovered 
with in-situ extraction methods. However, only 21.5 bil- 
lion barrels is considered economically recoverable with 
in-situ methods [15]. In-situ extraction is economically 
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viable when oil sands deposits are deeper than 100 m, 
and the deposits in Cold Lake are more than 400 m 
below the surface [14,16]. In-situ oil sands extraction op- 
erations in Alberta cover approximately 30 times the 
area that is available through surface mining, with more 
than twenty companies operating in the Cold Lake re- 
gion [14,15]. Cold Lake is also home to the Cold Lake 
First Nations, which includes three reserves, as well as 
the City of Cold Lake [15]. 

Much of the published research that examined PAH 
levels at Alberta oil sands extraction areas has been fo- 
cused on the effects of surface mining and bitumen up- 
grading in the Athabasca region. Consequently, the 
majority of this research has focused on regions with large 
open pit mining activities, and, in contrast, the environ- 
mental effects of in-situ extraction are relatively unknown. 
Kelly et al. [5,6] found increased heavy metals and PAH 
loadings in rivers downstream from bitumen upgrading fa- 
cilities in the Athabasca region, and concentrations of 
PAHs and heavy metals were found at elevated concentra- 
tions up to 50 km away from the upgrading facilities. 
More specifically, Kelly et al. [6] found an estimated 
annual loading of -1200 kg of particulate PAHs, 
with -500 kg of dissolved PAHs within a 25 km radius 
from the upgrading facilities. More recently, Kurek et al. 
[9] noted that PAH concentrations in lake sediments in the 
Athabasca region have increased concomitantly with com- 
mercial bitumen extraction operations at the Athabasca 
oil sands. The concentrations in all lakes have increased 
since development of the region began, with some dated 
lake sediment cores showing PAH concentrations increas- 
ing tenfold relative to predevelopment. Other studies in 
the Athabasca region have suggested that the region has 
high natural background concentrations that complicate 
differentiation between anthropogenic influences and 
those that would be expected as a result of natural pro- 
cesses such as erosion [2-4]. 

Many PAHs are known mutagens or carcinogens [17] 
that can enter the body through a number of routes of ex- 
posure including ingestion, inhalation, and dermal contact 
[18]. The larger 4-6 ringed PAHs, which are typically the 
most carcinogenic, are highly lipophilic and readily adsorb 
onto particulate matter in both water and aerosols, indi- 
cating that exposure to PAHs can occur via inhalation of 
atmospheric particulates, ingestion of contaminated par- 
ticulates (e.g., soil), or dermal contact with contaminated 
material (e.g., soil, sediment, water) [19-21]. Inadvertent 
ingestion of soil is frequently the dominant exposure 
pathway for soil at contaminated sites, particularly for 
non-volatile and semi-volatile contaminants such as 
carcinogenic PAHs. 

This study investigated the excess lifetime cancer risk 
for inhabitants of the Cold Lake region associated with 
exposures to PAHs via inhalation of airborne particulates 



and the inadvertent ingestion of soil. In conjunction with 
our companion quantitative soil ingestion study [22], the 
current work, which measured PAH contamination in 
both soil and air, provides a unique opportunity to quan- 
titatively evaluate the excess lifetime cancer risk to First 
Nations people and other residents in the Cold Lake re- 
gion. More specifically, measured soil ingestion rates for 
inhabitants that engage in wilderness activities [22], in 
conjunction with environmental PAH levels, are used to 
estimate the excess lifetime risk of cancer for residents 
in the Cold Lake region. Moreover, utilizing Human 
Health Risk Assessments (HHRA) methods currently ad- 
vocated by Health Canada [23], risk estimates for urban 
inhabitants of the region were also determined. We 
hypothesize that the excess cancer risk posed to inhabi- 
tants of the Cold Lake oil sands extraction region that 
are associated with potential exposures to carcinogenic 
PAHs in soil and air, constitutes an appreciable increase 
over background. 

Methods 

Study area 

The study was conducted near Cold Lake, Alberta, a site 
approximately 300 km northeast of Edmonton, and 
300 km southeast of the Athabasca oil fields. The Cold 
Lake region is one of the three major oil sands regions of 
Alberta, along with the Athabasca and Peace River de- 
posits. The region also contains the City of Cold Lake, 
with a population of approximately 14,000, the Canadian 
Forces Base Cold Lake, and four native reserves of the 
Cold Lake First Nations. The aboriginal people of the Cold 
Lake First Nations belong to the Dene Suline tribe, whose 
traditional lands ranged from south of Bonnyville to the 
northernmost point at Peter Pond Lake, Saskatchewan. 

Soil sampling and analysis 

This study collected 18 soil samples in August 2011 at 18 
locations in the Cold Lake region (Figure 1). Sampling 
involved collection of the surficial soil horizon with a pre- 
rinsed spatula, removal of surface vegetation, and place- 
ment into labeled WhirlPak™ bags. Soil samples were 
stored in a freezer at -20°C, and shipped to the laboratory 
in an ice packed cooler. The sampling locations are lo- 
cated within the areas where participants of the soil inges- 
tion study [22] spent significant portions of time. Soils 
were freeze-dried, lightly deconsolidated with a mortar 
and pestle, and sieved for 10 minutes to obtain the 63 \im 
fraction using an automated sieve shaker (Soil Test Engin- 
eering Model C1-592B, SoilTest, Evanston, IL, USA). The 
63 (im size fraction was chosen since this size fraction rep- 
resents the particulate matter size that adheres to hands, 
and is most likely to be inadvertently ingested [24] . 

To analyze soil samples for PAHs, 2 g of dried soil 
was mixed with Hydromatrix® (Agilent Technologies, 
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Figure 1 A map of the soil sampling sites in the Cold Lake region. Soil sampling sites are marked with a black square and numbered. Refer 
to Table 1. 



Mississauga, ON), and each sample was spiked with 
known concentrations of US EPA 16 PAH Cocktail 13 C 
(Cambridge Isotope Laboratories, Tewksbury, MA). 
PAHs were extracted from soil using accelerated solvent 
extraction (ASE 200 Accelerated Solvent Extraction Sys- 
tem, Dionex) at 100°C using 35% acetone: 65% hexane 
(Omnisolv grade solvents, Fisher Scientific, Ottawa, 
ON). Liquid-liquid extractions were used to separate 
polar extracts from non-polar fractions, which were col- 
lected using hexane, 2,2,4-trimethylpentane, and satu- 
rated Na 2 S0 4 . The non-polar fractions were purified to 
remove pigments and sulphur using USEPA method 
3640A with Envirogel columns on a preparative liquid 
chromatograph (Waters, Milford, Massachusetts; Agilent 
Technologies, Mississauga, ON). Purified samples were 
evaporated to 1 mL, and the PAH fraction separated on 
an activated 60-100 mesh Davisil 635 silica column 
(USEPA 3630C). The PAH fraction was analyzed with a 
6890 gas chromatograph and 5973 mass spectrometer 
by injecting 1 uL with pulsed splitless mode at 280°C on 
a DB-XLB 30 m x 0.18 urn x 180 urn column (Agilent 
Technologies, Mississauga, ON). The oven had an initial 
temperature of 60°C and was held for 2 minutes, and 
then increased at a rate of 6°C per minute up to 300°C, 
and held for 10 minutes. There was a constant flow of 
helium at a rate of 39 cm s" 1 for a total runtime of 52 mi- 
nutes. The mass spectrometer had a transfer line 



temperature of 280°C, a source temperature of 230°C 
and quadrupole temperature of 150°C. The 16 US EPA 
priority PAHs and 13 C labeled PAHs were analyzed and 
quantified with single ion monitoring. All PAH concen- 
trations were calculated using isotopic dilution with la- 
beled 13 C PAHs, method detection limits are provided in 
Additional file 1: Table SI. The 16 PAHs measured, and 
their potency equivalency factors (PEF) are provided in 
Additional file 2: Table S2. Potency Equivalency Factors 
(i.e., PEFs) were used to convert PAH levels to BaP 
equivalents for assessment of carcinogenic risk [25] . 

Air sampling and analysis 

A high volume continuous air sampler was setup on the 
Cold Lake First Nations reserve adjacent to the campsite 
used during the companion soil ingestion study [22]. The 
air sampler was calibrated with a 2100 magnehelic pres- 
sure gauge using a standard high-volume sampler calibra- 
tion procedure (Dwyer Instruments, Michigan City, IN). 
Flow through the high- volume sampler in m 3 min" 1 was 
calibrated against inches of H 2 0 magnehelic readings [26] . 
Glass fiber filters (GFF) (Whatman 110 mm 0.7 urn pore 
size, Fisher Scientific, Ottawa, ON) and large polyurethane 
foam (PUF) cartridges (Supelco, Oakville, ON) were re- 
placed every 24 hours. An average volume of 844.71 m 3 of 
air was sampled during each 24 hour period. Prior to use, 
all GFFs were placed in a muffle furnace for 6 hours at 
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500°C to remove organic carbon, and then pre-weighed 
and placed in sealed containers. Methanol rinsed metal 
tweezers were used to remove GFFs and PUF cartridges 
during sampling. After sampling, GFFS and PUF car- 
tridges were kept in a freezer at -20°C, and shipped to the 
laboratory in ice packed coolers. 

PAH analyses were conducted as described for soil sam- 
ples with the following differences. GFFs were weighed 
again, and placed in ASE cells that were filled with Hydro- 
matrix®, PAH method detection limits for GFF samples 
provided in Additional file 3: Table S3. PUF cartridges were 
cut in half, with the top and bottom halves analyzed separ- 
ately to determine breakthrough, and then placed in ASE 
cells that were filled with Hydromatrix®. PAH method de- 
tection limits shown in Additional file 4: Table S4. 

Risk assessment 

The risk assessment methodology advocated by both the 
US EPA and Health Canada was used for assessing excess 
lifetime cancer risk from exposure to PAHs [27,28]. Soil 
ingestion rates used to calculate the dose were determined 
through a First Nations soil ingestion study conducted 
over a 13 day period at Cold Lake [22]. Deterministic as- 
sessments to assess lifetime cancer risk through ingestion 
of soil were conducted using Eq. 1 [27], where Cs is the 
concentration in the soil (mg kg" 1 ), IR is the soil ingestion 
rate (mg d" 1 ), CF is the conversion factor (10" 6 ), EF is the 
exposure frequency, from 153 d yr" 1 to 365 d yr" 1 , ED is 
the exposure duration, which is 70 years to assess lifetime 
cancer risk, and BW is body weight (kg). Following the 
recommendations of the US EPA, 80 kg was used for this 
study [29]. AT is the averaging time of 25550 days (365 d 
y" 1 x 70 y) [27,30], and CSF is the cancer slope factor, 
which was 7.3 (mg kg" 1 d" 1 )" 1 for BaP ingestion [31]. Soil 
ingestion rates used to calculate dose were determined in 
the aforementioned 13-day soil ingestion rate study con- 
ducted at Cold Lake [22]. 



CsxIRxCF xEF x ED 

Risk = — — x CSF 

BW xAT 



(1) 



The assessment of excess lifetime cancer risk through 
inhalation of PAHs was conducted using Eq. 2 (from US 
EPA [32]). CA is the concentration in the air (mg m" 3 ), 
ET is the exposure time (24 hr d" 1 ), EF is exposure fre- 
quency (365 d yr" 1 ), ED is exposure duration (70 yr), and 
AT is an averaging time of 613200 hours (24 hr d" 1 x 
365 d y" r x 70 y), and IUR is the inhalation unit risk, 
which was 0.0033 (mg m" 3 )" 1 (from Health Canada [31]). 



CAxET xEF x ED 

Risk = x IUR 

AT 



(2) 



Monte Carlo simulation was employed to conduct 
probabilistic risk assessments following the US EPA [27] 



method. Eq. 1 was used with the same inputs as the deter- 
ministic risk assessments, except where indicated in the re- 
sults. Cs, IR, and EF were assumed to have lognormal or 
uniform distributions, and all Monte Carlo simulations 
were set for 10,000 iterations. Oracle™ Crystal Ball 
(11.1.2.2.0) software was used for Monte Carlo simulations. 

Excess cancer risk values were compared with literature 
values for four other petroleum extraction regions. More 
specifically, PAH concentrations for the 16 priority PAHs 
examined in this study were obtained from the litera- 
ture [33-35], converted to BaP equivalents using PEFs 
(Additional file 1: Table SI), and excess lifetime cancer risk 
assessed (Eq. 1) for maximum and mean BaP equivalent 
concentrations, assuming a 20 mg d" 1 soil ingestion rate. 

Health Canada defines excess lifetime cancer risk 
below 1E-05 as negligible, and this threshold was used 
to evaluate the results obtained in this study [28,31]. 

Results 

Soil PAH contamination 

Eighteen soil samples were collected and analyzed to de- 
termine concentrations of 16 priority PAHs and BaP 
equivalents at different distances from the oil sands facil- 
ities in Cold Lake, Alberta, and results are summarized in 
Table 1 and Figure 2. There was no significant relationship 
between PAH concentration, expressed as BaP equiva- 
lents, and distance from an oil sands facility (F (1> 16) = 1.73, 
p = 0.21). The maximum concentration of soil PAH was 
99.78 ng g" 1 of BaP equivalents sampled at 6.4 km 
(Figure 2, #14) from the nearest drilling pad. The second 
highest concentration of 79.00 ng g" 1 was observed at the 
roadside adjacent to one of the oil sands facilities. The 
mean soil concentration of BaP equivalents was 15.70 ± 
29.67 ng g" 1 , and the median concentration was 2.00 ng 
g" 1 . Most of the sampled soils had relatively low concen- 
trations with 13 soil samples below 5 ng g" 1 and 7 of these 
soil samples were below 1 ng g" 1 . The results obtained did 
not show a relationship between distance to in-situ oil 
sands facilities and soil PAH concentration. 

Atmospheric PAH contamination 

High-volume air sampling collected particulate and 
vapour-phase PAH samples over ten 24 hour sampling 
periods (Figure 3). The mean concentrations of BaP 
equivalents were 0.0043 ±0.013 ng m" 3 and 0.032 ± 
0.007 ng m" 3 for the particulate and vapour phases, re- 
spectively. The vapour phase PAH concentration, 
expressed as BaP equivalents, is an order of magnitude 
greater than the particle-associated PAH levels. There 
was no correlation between particulate and vapour- 
phase PAH concentrations. The maximum particulate 
PAH concentration of 0.042 ng m" 3 was obtained on the 
fifth day; the maximum gaseous PAH concentration of 
0.22 ng m" 3 was obtained on the tenth day. 
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Table 1 The 16 priority PAH concentrations, IPAH, and BaP equivalents, in soil samples collected in Cold Lake, Alberta 



Sample # 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


Distance from nearest pad (km) 


21.7 


0.1 


11.8 


0.1 


1.9 


3.3 


32 


30.6 


8.7 


34 


10.8 


18.9 


1.7 


6.4 


9.9 


26 


10.4 


0.1 


















PAH Concentration 


(ng g" 


) 














Naphthalene 


- 


1.9 


5.1 


- 


27.5 


- 


26.4 


- 


- 


- 


100.5 


- 


- 


- 


- 


- 


- 


- 


Acenaphthylene 


- 


- 


- 


38.2 


39.9 


30.4 


- 


24.7 


22.5 


22.3 


15.9 


22.4 


20.5 


6.3 


17.5 


30.9 


17.2 


93.5 


Acenaphthene 


- 


- 


- 


1.1 


1.9 


- 


- 


- 


0.5 


3.8 


- 


- 


18.1 


8.7 


- 


11 


39.7 


- 


Fluorene 


- 


0.6 


- 


- 


20.5 


- 


0.4 


- 


- 


- 


0.1 


1 


4.9 


2.7 


- 


- 


2.4 


1.4 


Phenanthrene 


- 


5.9 


- 


- 


23.5 


- 


1.4 


- 


- 


- 


1.9 


12.5 


29.2 


35.4 


4.3 


- 


33.1 


5.4 


Anthracene 


- 


- 


- 


- 


7.8 


- 


- 


- 


- 


- 


3.4 


- 


- 


- 


0.4 


2 


6.1 


0.6 


Fluoranthene 


- 


11.8 


- 


0.7 


7.6 


- 


4.2 


0.8 


- 


0.7 


13.8 


5.8 


5.9 


154.2 


18.4 


- 


14 


2 


Pyrene 


- 


46.5 


- 


- 


- 


- 


6.1 


2 


- 


- 


13.5 


2.6 


7 


66.1 


18.6 


- 


15.3 


6.2 


Benz[a]anthracene 


- 


4.4 


- 


- 


- 


- 


0.6 






0.6 


5.2 


0.6 


- 


5.3 


9 


- 


3.8 


0.1 


Chrysene 


0.7 


122.1 


0.3 


1 


4.3 


0.7 


5.2 


1.1 




2.7 


10.8 


6 


6 


33 


19.8 


- 


8.1 


1.2 


Benzo[b]fluoranthene 


0.6 


46.2 


- 


0.9 


- 


0.4 


1.7 






1 


10.5 


2.9 


3.5 


52.9 


25.3 


- 


8.6 


3.3 


Benzo[k]fluoranthene 




11.8 


















4.2 


0.3 




196.5 


9.5 




1.9 




Benzo[a] pyrene 




18.8 










0.6 






0.6 


7.1 


0.6 




17.9 


19.1 




3.8 




IndenoD ,2,3-cd]pyrene 




26.2 


















6.5 






66.8 


16.3 




3.3 




Dibenz[a,h]anthracene 




8 


























2.3 








Benzo[g,h,i]perylene 


1.9 


35.4 


















5.7 






20.9 


16.4 




4.8 


2.6 


IPAH 


3.2 


339.6 


5.4 


41.9 


133 


31.5 


46.6 


28.6 


23 


31.7 


199.1 


54.7 


95.1 


666.7 


176.9 


43.9 


162.1 


116.3 


BaP Equivalents 


0.7 


79.0 


3.0E-03 


1.3 


0.6 


0.7 


2.4 


0.3 


0.2 


1.5 


19.7 


3.5 


3.8 


99.8 


50.7 


0.3 


13.8 


4.3 



Detected concentrations in samples that were below the minimum detection limit had a value substituted with half the detection limit. Values that were not 
detected are assumed 0 and marked with a -. 
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Figure 2 The concentration of PAHs, expressed as BaP equivalents, in soil samples collected at Cold Lake, Alberta as a function of 
distance from the nearest oil pad. 
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Figure 3 Atmospheric concentrations of PAHs expressed as benzo[c/]pyrene equivalents. Panel A shows particulate PAH concentrations, 
and Panel B shows vapor-phase PAH concentrations. 



Cancer risk via soil ingestion 

Our earlier work [22] indicated that the 95th percentile soil 
ingestion rate is 361 mg d 1 , and the median rate is 37 mg 
d" 1 the former determined using Si as the soil tracer and the 
latter determined using the mean of Al and Si soil tracers. 
Using a deterministic risk assessment approach, the calcu- 
lated excess lifetime cancer risk at the maximum PAH con- 
centration observed in this study was 3.29E-06 at the 95th 
percentile soil ingestion rate, assuming exposure for 365 days 
per year [28]. In a Canadian climate, daily exposure for a full 
year seems unlikely because of the relatively harsh winter. A 
second exposure frequency of 153 days a year was also 
employed for the risk assessment. This value represents daily 
exposure for the period from May through September, and 
the excess lifetime cancer risk associated with this exposure 
at 361 mg d" 1 soil ingestion rate, and the maximum BaP 
equivalent soil concentration provided a value of 1.38E-06. 
Thus, there is no significant increase in the risk of cancer at 
this high soil ingestion rate; and moreover, the median soil 
ingestion rate the lifetime cancer risk was still substantially 
below 1E-05 (i.e., 1.41E-07). Furthermore, the mean BaP 
equivalent concentration of the collected Cold Lake soil 
samples (i.e., 15.7 ng g" 1 ) also yielded risk values that are 
substantially below 1E-05 (i.e., 5.17E-07 at the 95th soil in- 
gestion rate of 361 mg d" 1 and 365 day yr 1 exposure). 

Not surprisingly, the soil ingestion rate employed has a 
large effect on calculated excess lifetime cancer risk. At an 



assumed 365 days per year exposure (Figure 4) a signifi- 
cant increase in cancer risk occurs at a soil concentration 
of 304 ng g" 1 BaP equivalents for the 95th percentile soil 
ingestion rate. If lifetime cancer risk is calculated using the 
median soil ingestion rate of 37 mg d" 1 , lifetime cancer 
risk would meet the 1 case per 100,000 threshold at a BaP 
equivalent soil concentration of 2962 ng g" 1 . The adult soil 
ingestion rate of 20 mg d" 1 recommended by Health 
Canada for HHRA of contaminated sites, would require a 
soil concentration of 5479 ng g" 1 of BaP equivalents to 
reach an increased cancer risk of 1 case per 100,000 
people. Predictably, with an assumed exposure frequency 
of 153 days per year (Figure 5), which would represent 
non-winter months from May through September only, 
excess lifetime cancer risk in excess of 1E-05 would only 
be observed at much higher PAH concentrations. Using 
an EF of 153 day yr" 1 , at soil ingestion rates of 361 mg d" 1 , 
37 mg d" 1 , and 20 mg d" 1 , the 1.0E-05 cancer risk thresh- 
old occurs at BaP equivalent concentrations of 724 ng g" 1 , 
7066 ng g" 1 and 13072 ng g" 1 respectively. 

The BaP oral slope factor used in this study provides 
estimates of the risk of gastric (stomach) cancer through 
ingestion. In Alberta, Canada annual age-standardized 
new cases of stomach cancer are 7 per 100,000 [36]. 
Therefore, even at the maximum PAH concentration in 
this study, the median soil ingestion rate of 37 mg d" 1 
does not significantly increase the risk of gastric cancer 
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Figure 4 Calculated excess life time cancer risk as a function of soil PAH concentration (expressed as BaP equivalents) for different soil 
ingestion rates. Assumed exposure frequency of 365 days per year. The solid grey line denotes an excess lifetime risk of 1 extra cancer case per 
100,000 people. This 1E-05 risk occurs at a BaP concentration of 304 ng g" 1 , 2962 ng g" 1 , and 5479 ng g" 1 for soil ingestion rates of 361 mg d" 1 , 
37 mg d" 1 , and 20 mg d" 1 , respectively. The solid black vertical line denotes the maximum BaP equivalent concentration. Two of the chosen soil 
ingestion rates are the 95th percentile (361 mg d" 1 ) and median rates (37 mg d" 1 ) from the Cold Lake soil ingestion study. The third value 
(i.e., 20 mg d" 1 ) is the ingestion rate recommended by Health Canada for adult HHRA of contaminated sites. 



beyond the current background level (i.e., 7 cases per 
100,000). Even at higher soil ingestion rates, there is no 
significant increase in excess cancer risk (Table 2). Using 
the 90th percentile soil ingestion rate of 152 mg d" 1 and 
assumption of 365 d yr 1 exposure, the risk assessment 
predicts less than 0.1 new cancer cases per 100,000 over 
existing background. Using the 95th percentile soil inges- 
tion rate of 361 mg d" 1 , the predicted excess risk is still less 
than 0.1 additional case per 100,000 people for an EF of 
153 d yr" 1 . Thus, the results obtained, which are summa- 
rized in Table 2, indicate that the increased risk of gastric 
cancer associated with the median and mean PAH con- 
centrations found in the soils from the Cold Lake region is 
low and may be considered negligible. 

The excess cancer risk estimates shown above were 
compared to risk estimates based on published PAH 



levels for four other oil regions, using the Health Canada 
adult soil ingestion rate for HHRA of 20 mg d" 1 . The re- 
sults for the Niger Delta, Nigeria; Texas, USA; Vojvodina, 
Serbia; and Mathura, India and are summarized in 
Table 3. Although these regions had much higher PAH 
concentrations than the Cold Lake region, there calcu- 
lated excess cancer risk could still be termed negligible. 
The Texas, USA location had the greatest cancer risk 
level with 1.12E-05 for exposure at the maximum mea- 
sured BaP equivalent concentration, and the only loca- 
tion above the Health Canada 1E-05 threshold. It is 
important to note that the US EPA typically uses a 1E- 
06 cancer risk threshold, and using this threshold, the 
calculated excess lifetime cancer risk values for Vojvodina, 
Serbia, and Niger Delta, Nigeria, as well as the Texas, 
USA location. 
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Figure 5 Calculated excess life time cancer risk as a function of soil PAH concentration (expressed as BaP equivalents) for different soil 
ingestion rates. Assumed exposure frequency of 153 days per year. The solid grey line denotes an excess lifetime risk of 1 extra cancer case per 
100,000 people. This 1E-05 risk occurs at a BaP concentration of 724 ng g~ 1 , 7066 ng g~ 1 , and 13072 ng g" 1 for soil ingestion rates of 361 mg d" 1 , 
37 mg d" 1 , and 20 mg d" 1 , respectively. The solid black vertical line denotes the maximum BaP equivalent concentration. Two of the chosen soil 
ingestion rates are the 95th percentile (361 mg d" 1 ) and median rates (37 mg d" 1 ) from the Cold Lake soil ingestion study. The third value 
(i.e., 20 mg d" 1 ) is the ingestion rate recommended by Health Canada for adult HHRA of contaminated sites. 



Table 2 Rates of gastric cancer (per 100,000 people) in 
Alberta, and the calculated increases in cancer incidence 
associated with exposure to contaminated soils in Cold 
Lake, Alberta 

Stomach cancer New cases per 100,000 



Age-standardized incidence 3 


7 


Age-standardized mortality 3 


4 


Maximum PAH exposure 361 mg/d and 


7 


365 days/yr exposure frequency 




Maximum PAH exposure 361 mg/d and 


7 


153 days/yr exposure frequency 




Maximum PAH exposure 152 mg/d and 


7 


365 days/yr exposure frequency 




Maximum PAH exposure 152 mg/d and 


7 


153 days/yr exposure frequency 




Maximum PAH exposure 37 mg/d and 


7 


365 days/yr exposure frequency 





a: [32]. 

Calculated cases per 100,000 are the sum of the current background rates and 
the additional cases associated with exposure to contaminated soil. 



Probabilistic risk assessment 

Deterministic risk assessments (Figures 4 and 5) have 
been criticized for being too conservative (i.e., unneces- 
sarily inflate risk). To account for the uncertainty of in- 
put variables, probabilistic risk assessments have been 
recommended [27,38]. Excess lifetime cancer risk was 
assessed using the probabilistic Monte Carlo simulation 
approach (Figure 6), the BaP equivalent soil concentra- 
tions, and the soil ingestion rates determined in the Cold 
Lake soil ingestion study [22]. The calculated lifetime 
cancer risk in the Cold Lake region were determined to 
be negligible with 95% of the population associated with 
a risk level almost tenfold below the 1E-05 threshold (in- 
put variables in Table 4 Excess lifetime cancer was also 
calculated using the maximum sampled BaP equivalent 
soil concentrations (Figure 7), and the aforementioned 
soil ingestion rates from the Cold Lake study (input vari- 
ables Table 5). Again, calculated lifetime cancer risk is 
below the 1E-05 threshold, and thus negligible. A sensi- 
tivity analysis using Spearman's rank correlation coeffi- 
cient was carried out for the Monte Carlo risk analysis 
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Table 3 Comparison of BaP equivalent concentrations and lifetime cancer risk estimates for different locations 
reported in the literature 



Location 


n 


BaP Equivalent 3 (mg kg" 1 ) 


Lifetime cancer risk 


Maximum 


Mean 


Maximum 


Mean 


Niger Delta, Nigeria 6 


8 


5.3 


2.0 


4.05 E-06 


1 .50E-06 


Texas, USA C 


3 


14.7 


11.0 


1.12E-05 


8.41 E-06 


Vojvodina, Serbia d 


7 


8.2 


2.6 


6.28E-06 


1 .98E-06 


Mathura, India 6 


29 




0.34 




2.59E-07 


Cold Lake, Canada 


18 


0.1 


0.02 


7.63E-08 


1 .20E-08 



a: PAH concentrations converted to BaP equivalent concentrations using toxic equivalency factors [24]. 
b: [34]. 
c: [33]. 
d: [35]. 
e: [37]. 

Cancer risk was assessed using the maximum Cs, an ED of 70 years, an EF of 153 days, BW of 80 kg, an AT of 25550 days, and CSF ora | of 7.3 (mg kg" 1 d" 1 )" 1 . 



equation used for the probabilistic risk analysis [38]. The 
three variables examined were Cs, EF, and IR (Table 6). 
All three variables were positively correlated with life- 
time cancer risk, with ingestion rate having the largest 
impact on lifetime cancer risk, closely followed by Cs, 
with EF having the smallest influence on risk. 



Cancer risk via inhalation 

Atmospheric samples of PAHs were collected near the 
base camp, approximately 10 km from the in-situ mining 
facilities, and inhalation rates recommended by Health 
Canada were used to assess excess lifetime cancer risk 
(Table 7) to recorded daily concentrations of BaP 




mean = 0.314- 95% = 1 .31 



Risk (10-7) 

Figure 6 Cumulative probability distribution of a Monte Carlo simulation based on 10,000 trials and the input variables from Table 4. 

The soil concentrations are the distribution of BaP equivalents measured at Cold Lake (Figure 1), the IR values were from Cold Lake soil ingestion 
study [22], and the exposure frequency is 153 days. The mean risk level is 3.42E-8, and the 9 % risk level is 1.31E-7. 
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Table 4 Input variables for Monte Carlo simulation 
displayed in Figure 6 



Input variable 


Distribution 


Parameters 


Cs (mg/kg d.w. BaP 
equivalent) 


Lognormal 


Mean = 0.01 57, SD = 0.0297 


IR (mg/d) 


Lognormal 


Mean = 52, SD= 119 


EF d/yr 


Uniform 


Min = max= 153 


ED yr 


Uniform 


Min = max = 70 


BW kg 


Uniform 


Min = max = 80 


AT 


Uniform 


Min = max = 25550 


CSF ora | (mg/kg-day) 


Uniform 


Min = max = 7.3 



equivalents. The calculated excess lifetime cancer risk 
was far below the 1E-05 threshold, even at the max- 
imum measured atmospheric BaP equivalent concentra- 
tion of 2.15E-04 ng m" 3 excess lifetime cancer risk 
through inhalation was 7. IE- 10. 

Discussion 

Soil and atmospheric PAH contamination 

This study was designed to assess the excess lifetime 
cancer risk posed by exposures to PAHs in a region of 
Alberta that contains the worlds largest in-situ bitumen 



extraction operation, a process expected to overtake 
other conventional forms of petrochemical production. 
It was expected that PAH concentrations would follow a 
similar trend to that seen in Kelly et al. [6], where PAH 
concentrations were elevated up to 50 km from the oil 
sands upgraders in the Athabasca oil sands region. Al- 
though there is in-situ extraction in the Athabasca area, 
the region is dominated by bitumen extraction through 
open pit mining. Furthermore, the several upgrading and 
refining facilities in the Athabasca region are suspected 
of being major sources of atmospheric PAHs [6]. The 
Cold Lake region, despite being a large oil sands region, 
has a very different industrial landscape compared to the 
Athabasca region. The oil sands operations extract all 
bitumen using in-situ steam extraction methods, and 
the extracted bitumen is transported via pipeline to 
Lloydminster or Edmonton for upgrading and refining. 
The focus on in-situ extraction, and the lack of upgrading 
or oil refining in the Cold Lake region, may account for 
the comparatively low PAH concentrations observed in 
Cold Lake. The soil PAH concentrations recorded in this 
study suggest that in-situ oil sands development in the 
Cold Lake region do not make substantial contributions 
to local increases in soil PAH concentrations. 
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Figure 7 Cumulative probability distribution of a Monte Carlo simulation based on 10,000 trials and the input variables from Table 5. 

The soil concentration was the highest BaP equivalent concentration measured at Cold Lake (Figure 1), IR values were from Cold Lake soil 
ingestion study [22], and the exposure frequency is 153 days. The mean risk level is 1.82E-7, and the 95% risk level is 6.72E-7. 
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Table 5 Input variables for Monte Carlo simulation 



displayed in Figure 7 



Input variable 


Distribution 


Parameters 


Cs (mg/kg d.w. BaP equivalent) 


Uniform 


Min = max = 0.0998 


IR (mg/d) 


Lognormal 


Mean = 52 SD= 119 


EF d/yr 


Uniform 


Min = max= 153 


ED yr 


Uniform 


Min = max = 70 


BW kg 


Uniform 


Min = max = 80 


AT 


Uniform 


Min = max = 25550 


CSF ora | (mg/kg-day) 


Uniform 


Min = max = 7.3 



Risk assessment 

The results obtained failed to show an excess lifetime 
risk of cancer above the 1E-05 threshold, using the 95th 
percentile soil ingestion rate of 361 mg d" 1 [22] (Figure 4 
and Table 2) and a 365 day per year EF. This EF will al- 
most certainly contribute to an unrealistically inflated 
cancer risk due to the likely changes in seasonal soil in- 
gestion rates, (i.e., an individual with an IR SO ii at the 95th 
percentile would be expected to have a much lower IR SO ii 
during winter months). Unfortunately, no published soil 
ingestion rate studies have investigated seasonal varia- 
tions in ingestion rates. Thus, the likelihood, magnitude 
and direction of a seasonal change in IR SO ii is unknown. 
It is reasonable to assume that soil ingestion rates during 
warmer months, when First Nations' people who follow 
a traditional lifestyle spend greater amounts of time out- 
side, would be substantially higher than during the win- 
ter months. Nonetheless, the estimated cancer risk 
associated with the reduced EF of 153 days per year is 
still below the 1E-05 threshold. 

The Canadian Council of Ministers of the Environment 
(CCME) has set soil quality guidelines for BaP at 700 ng 
g 1 for the protection of human health [39,40]. The guide- 
lines are based on an adult soil ingestion rate of 20 mg d 1 , 
and this ingestion rate may not be applicable to all popula- 
tions. Although the 37 mg d" 1 median soil ingestion rate 
from the Cold Lake First Nations soil ingestion study is al- 
most twice the recommended IR soi i for adults, the risk es- 
timates presented here indicate that the excess lifetime 
cancer risk associated with levels below the CCME soil 
quality guidelines and the higher IR soi i are below 1E-05, 
and may be regarded as negligible. Probabilistic risk as- 
sessments give a better understanding of the cancer risk 



Table 6 Sensitivity analysis Spearman's rank correlation 
coefficient values 



Input variable 


Coefficient 


Cs mg kg" 1 


0.61 


EF d yr" 1 


0.11 


IR mg d" 1 


0.71 



Table 7 Cancer risk estimates associated with inhalation 
of atmospheric PAHs in Cold Lake, Alberta 



BaP equivalent PAH concentration ug m" 3 Lifetime cancer risk 



2.15E-04 (maximum) 


7.10E-10 


1.32E-08 (minimum) 


4.37E-14 


3.68E-05 (mean) 


1.20E-10 


7.08E-06 (median) 


2.34E-1 1 


10~ 6 . 




10" 7 . 





posed to a population, as opposed to simple point esti- 
mates based on measured Cold Lake soil ingestion rates 
(Figures 6 and 7). Nevertheless, probabilistic assessments 
also indicated that the excess lifetime cancer risk associ- 
ated with exposures to soil PAH levels observed at Cold 
Lake are well below 1E-05, and can be regarded as negli- 
gible. The significance of this determination of negligible 
risk should not be overlooked since in-situ bitumen ex- 
traction is expected to overtake open pit mining as the 
dominant form of bitumen extraction in Alberta [15]. 

The calculated increases in cancer incidence that would 
be expected from inadvertent ingestion of contaminated 
soil at Cold Lake, even at the 95% IR soi i rate and maximum 
BaP equivalent soil concentrations (Table 2), are well below 
the background incidence of gastric cancer (i.e., 7E-05), 
and thus, the expected increase over background can also 
be considered to be negligible [36]. In contrast, a study that 
examined excess lifetime cancer risk from oral exposures 
to PAHs in settled house dust (i.e., inadvertent ingestion) 
noted far higher risk estimates. More specifically, for the 
maximum PAH concentration recorded in the study by 
Maertens et al. [41], and a 50 mg d" 1 dust ingestion rate, 
the authors noted an increase of 27.4 extra gastric cancer 
cases per 100,000, compared to less than 1 excess per 
100,000 people in this study for the maximum PAH con- 
centration. Similarly, Williams et al. [42] investigated the 
excess cancer risk from exposure to settled house dust and 
soil at locations adjacent to parking lots sealed with coal- 
tar containing products. Using an IR distribution with a 
mean of 27 mg d" 1 , exposure to PAHs in settled house dust 
yielded an increase of 1.8 cases per 100,000 people, and 1.2 
cases per 10,000 at the 50th and 95th percentiles, respect- 
ively. Exposure to soils adjacent to parking lots sealed with 
coal-tar containing products were associated with an even 
a greater increase of 7.3 cases per 100,000 at the 50th per- 
centile, and 4.3 cases per 10,000 people at the 95th per- 
centile. Although the studies by Maertens et al. [41] and 
Williams et al. [42] examined very different environments, 
they reinforce the notion that the calculated risks attribut- 
able to PAHs in soils near an in-situ oil sands facility might 
be reasonably regarded are indeed negligible. 

It should be noted that other studies of excess lifetime 
cancer risk attributable to ingestion of PAHs in contami- 
nated soil and street dust have also yielded values that 
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can be considered negligible. For example, a study by 
Wang et al [42] that calculated risk from exposures to 
PAHs in urban street dust revealed a mean risk from 
dust ingestion of 2.51E-6. The mean BaP equivalent con- 
centration from the study by Wang et al [42] was 
0.47 mg kg" 1 , which is much greater than the average for 
Cold Lake soils of (i.e., 0.016 mg kg" 1 ). Man et al. [43] 
used PAH levels in 12 different soils collected in Hong 
Kong to assess the excess lifetime cancer risk. At the 
median level, no soil yielded a cancer risk above the IE- 
05 threshold, and only one site, near a car dismantling 
workshop yielded an excess cancer risk at the 95th per- 
centile of 24 per 100,000. Among the other sites in the 
Man et al. [43] study that might be expected to be asso- 
ciated with cancer risk values greater than the 1E-05 
threshold was an open burning site where large amounts 
of pyrogenic PAHs would be expected. Nevertheless, the 
calculated rates were below the 1E-05 threshold. Al- 
though only one location from the Man et al. study 
yielded cancer risk estimates above the 1E-05 threshold, 
PAH levels at seven of the sampled sites, which ranged 
from 0.090 mg kg" 1 BaP equivalents to a maximum of 
4.93 mg kg" 1 near a car dismantling workshop, were 
above the Cold Lake average of 0.016 mg kg" 1 . 

Despite the fact that the deterministic scenarios 
yielded risk values below the 1E-05 threshold, the deter- 
ministic risk assessments may still be expected to inflate 
risk values over the true lifetime cancer risk. In contrast, 
the probability distribution functions (PDF) (Figure 6) 
are likely to be far more representative of lifetime cancer 
risk in the Cold Lake region, since the PDF functions 
can account for the distribution of some or all the vari- 
ables included in a risk assessment calculation [38]. For 
example, the soil ingestion rate likely varies across a 
population, as evidenced by the relatively high standard 
deviations noted in the companion soil ingestion study 
[22]. Additionally, PAH concentrations are also highly 
variable (Figure 2 and Table 1). The same risk assess- 
ment was run with both a deterministic and probabilistic 
approach. Although the deterministic assessment, which 
used an IR SO ii of 37 mg d" 1 and the mean soil PAH con- 
centration, yielded a cancer risk of 2.22E-8, the probabil- 
istic approach, which used a distribution for PAH soil 
concentration and IR SO ii (Table 4), yielded a very similar 
mean cancer risk of 3.37E-8. Nevertheless, it is critical to 
note that the probabilistic assessment and the cumula- 
tive probability figure (i.e., Figure 6) provide an ability to 
determine that the mean risk level of 3.39E-8 accounts 
for more than 80% of the population. Thus, the expected 
probabilistic cancer risk level for most of the population 
is even lower than the level calculated using the deter- 
ministic equation [38]. In the absence of a PDF, the cal- 
culated risk may only represent a small part of the 
population. 



Few published studies report soil PAH concentrations 
in petroleum extraction and refining regions; however, 
the concentrations that have been recorded are much 
greater than those found in the Cold Lake oil sands re- 
gion (Table 3). More specifically, the mean BaP equiva- 
lent concentrations observed at other sites are at least 
an order of magnitude greater than those observed in 
this study, and only the mean concentration from a site 
in India, which reported a BaP equivalent concentration 
of 0.34 mg kg" 1 , was below the Canadian CCME guide- 
line value [37] . The cancer risk estimates associated with 
the PAHs recorded at other petroleum extraction and 
refining locations are also much greater than those cal- 
culated for Cold Lake, but nevertheless, still below the 
Health Canada threshold. However, it should be empha- 
sized that the industrial landscape at Cold Lake is different 
from the other examined locations. For example, with the 
exception of the Texas location, where soils were contami- 
nated by faulty storage facilities, all the other locations 
with high PAH levels were in close proximity to refinery 
operations. This suggests that atmospheric release and de- 
position in the vicinity of oil refining operations contrib- 
utes to regional soil contamination. 

The sensitivity analysis (Table 6) shows that the soil in- 
gestion rate metric has the largest effect on cancer risk es- 
timation [27]. Sensitivity analysis data in other published 
studies are rare, but one study that also assessed cancer 
risk from PAHs through inhalation and dermal exposure, 
also found EF to have the lowest influence on lifetime can- 
cer risk, relative to other variables [44]. Williams et al. 
[45] also conducted a sensitivity analysis on the Cs and IR 
(dust and soil) variables, and found that, for adults, the Cs 
parameter had the greatest influence on assessments that 
examined exposures to PAHs in soil or house dust adja- 
cent to pavement where coal-tar containing products were 
not used. In contrast, if coal-tar containing sealant was 
used, then IR soi i had a larger influence on cancer risk esti- 
mation. These results are consistent with this study, where 
both IR and Cs were found to have similar influence on 
cancer risk estimation. 

The cancer slope factor employed in this study is rec- 
ommended by Health Canada. However, it is important 
to note that other jurisdictions and regulatory agencies 
may use different slope factors for determinations of 
excess lifetime cancer risk. For example, the California 
Office of Environmental Health Hazard Assessment 
(OEHHA) uses an oral slope factors of 2.9 (mg kg" 1 d" 1 )" 1 , 
which is markedly different from the Health Canada 
value of 7.3 (mg kg" 1 d" 1 ). Nevertheless, with respect to 
the risk assessments conducted in this study, modest dif- 
ferences in slope factors would not significantly alter the 
excess lifetime cancer risk values associated with PAH 
exposures in Cold Lake (i.e., risk values still well below 
the 1E-05 threshold). 
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Cancer risk via inhalation 

The atmospheric concentrations of PAHs are also asso- 
ciated with risk estimates that are well below the thresh- 
old value of 1E-05. One published study reported rural 
atmospheric BaP equivalent concentrations in Canada as 
a median of 0.14 ng m" 3 , with major cities such as 
Toronto and Winnipeg having median BaP equivalent 
concentrations of 076 ng m~ 3 and 0.14 ng m" 3 , respectively, 
compared to the maximum BaP equivalent concentra- 
tion measured in this study of 0.21 ng m" 3 [20]. In com- 
parison, the World Health Organization recommends an 
atmospheric BaP concentrations of 0.12 ng m" 3 to ensure 
that excess risk does not exceed 1 cancer case per 
100,000 people [19,46], current guidelines from the 
Alberta Government require an annual average BaP atmos- 
pheric concentration of 0.30 ng m" 3 [47], and the 
Ontario Ministry of the Environment (MOE) specifies 
an Ambient Air Quality guideline of 1.1 ng m" 3 over a 
24 hour period [48]. Atmospheric BaP concentrations 
measured in the Cold Lake region are also typically 
lower than levels in major urban centers of Alberta. 
WBK & Associates Inc. [48] reported a 7 year median 
concentration of 0.12 ng m" 3 for Calgary and 0.14 ng 
in" 3 for Edmonton, with maximum annual average concen- 
trations of 0.37 ng m" 3 and 0.41 ng m" 3 in Calgary and 
Edmonton, respectively. Only one of the observations 
recorded in the current is similar to typical urban BaP 
concentrations in Alberta, and the atmospheric concen- 
tration in Cold Lake can be regarded as fairly represen- 
tative of a rural region in Alberta [18,49]. 

The low atmospheric PAH concentrations yielded life- 
time cancer risk estimates from inhalation of atmos- 
pheric PAHs that are negligible for all concentrations, 
and even below the more conservative US EPA cancer 
risk threshold of 1E-6 (Table 7). It is difficult to know if 
the atmospheric PAH levels recorded over the 10 day 
sampling period are representative; nevertheless, the cal- 
culated risk estimates are negligible, and other studies 
that compared excess cancer risk from PAH exposures 
via inhalation and ingestion have noted that inhalation 
exposures make relatively small contributions to total 
risk [42,49]. 

Conclusions 

This study examined excess lifetime cancer risks posed 
by exposures to PAHs in soil and air at an area close to 
the worlds largest in-situ bitumen extraction area. The 
deterministic and probabilistic risk assessment, which 
employed IR SO ii values from the companion First Nations 
soil ingestion study [22], both indicate that cancer risk is 
below the 1E-05 threshold, and can thus be considered 
negligible in most jurisdictions. Moreover, the PDF 
showed that 95% of the population is below the risk 
threshold of 1E-05. Additionally, the risk associated with 



PAH exposures via inhalation was also considered negli- 
gible. The hypothesis that there is an increased excess 
lifetime cancer risk from PAH exposures posed to First 
Nations people in the Cold Lake region is not supported 
by this study. 

Additional files 



Abbreviations 

ASE: Accelerated solvent extraction; AT: Averaging time; BaP: Benzo[a]pyrene; 
BW: Body weight; CCME: Canadian Council of Ministers of the Environment; 
CF: Conversion factor; CS: Concentration; CSF: Cancer slope factor; 
ED: Exposure duration; EF: Exposure frequency; GFF: Glass fibre filter; 
HHRA: Human health risk assessment; IR: Intake rate; PAH: Polycyclic aromatic 
hydrocarbon; PEF: Potency equivalency factors; PDF: Probability density 
functions; PUF: Polyurethane foam cartridge; SD: Standard deviation. 

Competing interests 

The authors have no conflict of interest to declare. 
Authors' contributions 

Gl carried out the collection of soil and atmospheric samples, processed all 
samples, performed the data analyses, and prepared a draft of the 
manuscript. JMB conceived of the study, participated in its design, and 
helped draft the manuscript. JRD participated in study design, data analysis, 
and helped draft the manuscript. LEK participated in study design, method 
development and data analysis. PAW participated in study design, data 
analysis, and helped draft the manuscript. All authors have read and 
approved the final manuscript. 

Acknowledgements 

We would like to thank Cold Lake First Nation for the use of their land 
during the study. As well as Tribal Chiefs Ventures Inc. for their support and 
help in coordinating the field work. Thank you to Elliott Skierszkan for 
laboratory assistance. Funding for this research was provided by Health 
Canada under the First Nations Environmental Contaminants Program. 

Received: 19 July 2013 Accepted: 7 February 2014 
Published: 12 February 2014 

References 

1. CIA: CIA - The world factbook. https://www.cia.gov/library/publications/the- 
world-factbook/rankorder/2244rank.html. 

2. Headley JV, Akre C, Conly FM, Peru KM, Dickson LC: Preliminary 
characterization and source assessment of PAHs in tributary sediments 
of the Athabasca River, Canada. Environ Forensics 2001, 2:335-345. 

3. Headley JV, Marsh P, Akre C, Peru KM, Lesack L: Origin of Polycyclic 
Aromatic Hydrocarbons in Lake Sediments of the Mackenzie Delta. 

J Environ Sci Health A Tox Hazard Subst Environ Eng 2002, 37:1 1 59-1 180. 



Additional file 1: Table SI. PAH method detection limit and 
calculation for soil samples. Ce is analyte concentration, Ve, is injection 
analyte volume, Vs is average sample mass measured, R% is the average 
recovery rate of 13 C labeled PAHs, and MDL is the calculated method 
detection limit. 

Additional file 2: Table S2. List of the 16 priority PAHs measured for 
this study and their Potency Equivalency Factors [24], 

Additional file 3: Table S3. Method detection limit and calculation for 
particulate PAH samples measured on GFFs. Ce is analyte concentration, 
Ve, is injection analyte volume, Vs is average sample mass measured, R% 
is the average recovery rate of 13 C labeled PAHs, and MDL is the 
calculated method detection limit. 

Additional file 4: Table S4. Method detection limit and calculation for 
gaseous PAH samples measured on PUF cartridges. Ce is analyte 
concentration, Ve, is injection analyte volume, Vs is average sample mass 
measured, R% is the average recovery rate of 13 C labeled PAHs, and MDL 
is the calculated method detection limit. 



Irvine et al. Environmental Health 2014, 13:7 
http://www.ehjournal.net/content/1 3/1/7 



Page 14 of 15 



4. Akre C, Headley JV, Conly FM, Peru KM, Dickson LC: Spatial patterns of 
natural polycyclic aromatic hydrocarbons in sediment in the lower 
Athabasca River. J Environ Sci Health A Tox Hazard Subst Environ Eng 2004, 
39:1163-1176. 

5. Kelly EN, Schindler DW, Hodson PV, Short JW, Radmanovich R, Nielsen CC: 
Oil sands development contributes elements toxic at low concentrations 
to the Athabasca River and its tributaries. Proc Natl Acad Sci USA 2010, 
107:16178-16183. 

6. Kelly EN, Short JW, Schindler DW, Hodson PV, Ma M, Kwan AK, Fortin BL: Oil 
sands development contributes polycyclic aromatic compounds to the 
Athabasca River and its tributaries. Proc Natl Acad Sci USA 2009, 
106:22346-22351. 

7. Wiklund JA, Hall Rl, Wolfe BB, Edwards TWD, Farwell AJ, Dixon G: Has Alberta 
oil sands development increased far-field delivery of airborne contaminants 
to the Peace-Athabasca Delta? Sci Total Environ 201 2, 433:379-382. 

8. Hall Rl, Wolfe BB, Wiklund JA, Edwards TWD, Farwell AJ, Dixon G: Has 
alberta oil sands development altered delivery of polycyclic aromatic 
compounds to the peace-athabasca delta? PLoS One 2012, 7:1-17. 

9. Kurek J, Kirk JL, Muir DCG, Wang X, Evans MS, Smol JP: Legacy of a half 
century of Athabasca oil sands development recorded by lake 
ecosystems. Proc Natl Acad Sci USA 201 3, 1 1 0:1 -6. 

10. Timoney KP, Lee P: Polycyclic aromatic hydrocarbons increase in 
Athabasca River Delta sediment: temporal trends and environmental 
correlates. Environ Sci Technol 201 1, 45:4278-4284. 

11. Timoney KP, Lee P: Does the Alberta tar sands industry pollute? The 
scientific evidence. Open Conservat Biol J 2009, 3:65-81 . 

12. CBC News: Cancer rates downstream from oilsands to be probed, http:// 
www.cbc.ca/news/canada/edmonton/cancer-rates-downstream-from- 
oilsands-to-be-probed-1 .1 030670. 

1 3. Weber B: Proposed oil sands health study derailed after aboriginal band 
pulls support over cancer rate debate, http://business.financialpost.com/ 
2013/02/22/oil-sands-health-survey-collapses-after-aboriginal-band-pulls- 
support/? Isa=67a5-c50c. 

14. Moorhouse J, Huot M, Dyer S: Drilling deeper: the in situ oil sands report card. 
AB: The Pembina Institute, Drayton Valley; 2010. 

15. Teare M, Burrowes A, Baturin-Pollock C, Rokosh D, Evans C, Gigantelli P, 
Marsh R: ST98-2012: Alberta's energy reserves 201 1 and supply/demand outlook 
2012-2021. Calgary, Alberta: Energy Resources Conservation Board; 2012. 

16. Gosselin P, Hrudey S, Anne Naeth M, Plourde A, Therrien R, Van Der Kraak G, 
Xu Z: Environmental and Health Impacts of Canada's Oil Sands Industry. 
Ottawa, ON: The Royal Society of Canada; 2010. 

1 7. Arey J, Beland FA, Borlak J, Burchiel SW, Cavalieri E, Feron VJ: IARC 
monographs on the evaluation of carcinogenic risks to humans Volume 92 
some non-heterocyclic polycyclic aromatic hydrocarbons and some related ex- 
posures. Lyon, France: International Agency for Research on Cancer; 2010. 

18. Stowers SJ, Anderson MW: Formation and persistence of benzo(a)pyrene 
metabolite-DNA adducts. Environ Health Perspect 1985, 62:31-39. 

19. Bostrom C, Gerde P, Hanberg A, Jernstrom B, Johansson C, Kyrklund T, 
Rannug A, Tornqvist M, Victorin K, Westerholm R: Cancer risk 
assessment, indicators, and guidelines for polycyclic aromatic 
hydrocarbons in the ambient air. Environ Health Perspect 2002, 
110(Suppl 3):451-488. 

20. Government of Canada, Environment Canada, Health Canada: Canadian 
Environmental Protection Act Priority Substances List Assessment Report - 
Polycyclic Aromatic Hydrocarbons. Ottawa, ON: Government of Canada, 
Environment Canada, Health Canada; 1994. 

21. Lima A, Farrington J, Reddy CM: Combustion-derived polycyclic aromatic 
hydrocarbons in the environment — a review. Environ Forensics 2005, 
6:109-131. 

22. Irvine G, Doyle JR, White PA, Blais JM: Soil ingestion rate determination in 
a rural population of Alberta, Canada practicing a wilderness lifestyle. Sci 

Total Environ 2014, 470-471:138-146. 

23. Wilson Scientific Consulting Inc, Meridian Environmental Inc: Critical Review 
of Soil Ingestion Rates for use in Contaminated Site Human Health Risk 
Assessments in Canada. Calgary, AB: Wilson Scientific Consulting Inc, 
Vancouver, BC & Meridian Environmental Inc.; 2006. 

24. Doyle JR, Blais JM, White PA: Mass balance soil ingestion estimating 
methods and their application to inhabitants of rural and wilderness 
areas: a critical review. Sci Total Environ 2010, 408:2181-2188. 

25. Schneider K, Roller M: Cancer risk assessment for oral exposure to PAH 
mixtures. J Appl Toxicol 2002, 22:73-83. 



26. St-Amand AD, Mayer PM, Blais JM: Modeling atmospheric vegetation 
uptake of PBDEs using field measurements. Environ Sci Technol 2007, 
41:4234-4239. 

27. Bennet DA, Chang SS, Cooper DE, Dinan J, Hofmann EL, Barry TM, Burke A, 
Galizia A, Olsen M, Jafolla NR, Simon TW, Thorns SR, Mucha A, Chapman J, 
Martinez ML, Griffin S, Henningsen G, Hoff D, Goulet J: Risk Assessment 
Guidance for Superfund: Volume 3 - Part A, Process for Conducting Probabilistic 
Risk Assessment. Washington: U,S. Environmental Protection Agency; 2001. 

28. Health Canada: Federal contaminated site risk assessment in Canada - Part I: 
Guidance on human health preliminary quantitative risk assessment (PQRA) 
Version 2.0. Ottawa, ON: Health Canada; 2012. 

29. Moya J, Phillips L, Schuda L, Wood P, Diaz A, Lee R, Clickner R, Birch RJ, 
Adjei N, Blood P, Chapman K, de Castro R, Mahaffey K: Exposure Factors 
Handbook: 201 1 edition. U.S. Environmental Protection Agency: Washington; 
2011. 

30. Kumar B, Gaur R, Kumar S, Sharma CS: Environmental and human health risk 
assessment of Benzo(a)pyrene levels in agricultural soils from the National 
Capital Region, Delhi, India. Hum Ecol Risk Assess 2013, 19:1 18-125. 

31 . Health Canada: Federal contaminated site risk assessment in Canada - Part II: 
Health Canada toxicological reference values (TRVs) and chemical-specific fac- 
tors, Version 2.0. Ottawa, ON: Health Canada; 2010. 

32. Bailey M, Benson B, Crawford D, Draugelis A, Foos B, Foureman G, Griffin S, 
Hodoh O, Hubbard J, Johnson A, Johnson J, Koporec K, Levinson S, Murphy 
D, Schuver H, Sivak M, Stanek J, Stralka D, Taylor T, Whalan J, Winchester E, 
Overstreet C, Stiber N: Risk Assessment Guidance for Superfund: Volume I: 
Human Health Evaluation Manual (Part F, Supplmenetal Guidance for 
Inhalation Risk Assessment). Washington: U,S. Environmental Protection 
Agency; 2009. 

33. Bojes HK, Pope PG: Characterization of EPA's 16 priority pollutant 
polycyclic aromatic hydrocarbons (PAHs) in tank bottom solids and 
associated contaminated soils at oil exploration and production sites in 
Texas. Regul Toxicol Pharmacol 2007, 47:288-295. 

34. Duke O, Albert I: Spatial variation and distribution of polycyclic aromatic 
hydrocarbons in soil. B Chem Soc Ethiopia 2007, 21:331-340. 

35. Skrbic B, Durisic-Mladenovic N: Levels of PAHs in soil samples from the 
vicinity of oil refinery Novi Sad-Serbia. Kuwait J Sci Eng 2009, 36:63-75. 

36. Canadian Cancer Society's Steering Committee on Cancer Statistics: 
Canadian Cancer Statistics 2012. Ottawa, ON: Canadian Cancer Society, 
Toronto & Statistics Canada; 2012. 

37. Tiwari JN, Chaturvedi P, Ansari NG, Patel DK, Jain S, Murthy R: 
Assessment of Polycyclic Aromatic Hydrocarbons (PAH) and heavy 
metals in the vicinity of an oil refinery in India. So/7 Sediment Contam 
2011, 20:315-328. 

38. Oberg T, Bergback B: A review of probabilistic risk assessment of 
contaminated land. J Soil Sediment 2005, 5:213-224. 

39. Canadian Council of Ministers of the Environment: Canadian Soil Quality 
Guidelines for Carcinogenic and Other Polycyclic Aromatic Hydrocarbons 
(Environmental and Human Health Effects). Gatineau, QC: Canadian Council 
of Ministers of the Environment; 2010. 

40. Canadian Council of Ministers of the Environment: A Protocol for the 
Derivation of Environmental and Human Health Soil Quality Guidelines. 
Winnipeg, MN: Canadian Council of Ministers of the Environment; 2006. 

41. Maertens R, Yang X, Zhu J, Gagne R, Douglas G, White PA: Mutagenic and 
carcinogenic hazards of settled house dust I: Polycyclic aromatic 
hydrocarbon content and excess lifetime cancer risk from preschool 
exposure. Environ Sci Technol 2008, 42:1 747-1 753. 

42. Wang W, Huang M, Kang Y, Wang H, Leung AOW, Cheung KC, Wong MH: 
Polycyclic aromatic hydrocarbons (PAHs) in urban surface dust of 
Guangzhou, China: status, sources and human health risk assessment. 
Sci Total Environ 201 1 , 409:451 9-4527. 

43. Man YB, Kang Y, Wang HS, Lau W, Li H, Sun XL, Giesy JP, Chow KL, Wong 
MH: Cancer risk assessments of Hong Kong soils contaminated by 
polycyclic aromatic hydrocarbons. J Hazard Mater 2013:1-7. In Press. 

44. Chen S-C, Liao C-M: Health risk assessment on human exposed to envir- 
onmental polycyclic aromatic hydrocarbons pollution sources. Sci Total 
Environ 2006, 366:112-123. 

45. Williams ES, Mahler BJ, Van Metre PC: Cancer risk from incidental ingestion 
exposures to PAHs associated with coal-tar-sealed pavement. Environ Sci 
Technol 2013, 47:1101-1109. 

46. World Health Organization: Air quality guidelines for Europe. Secondth 
edition. Copenhagen: World Health Organization; 2000. 



Irvine et al. Environmental Health 2014, 13:7 
http://www.ehjournal.net/content/1 3/1/7 



Page 15 of 15 



47. Government of Alberta: Alberto Ambient Air Quality Objectives and Guidelines 
Summary. Edmonton, AB; 2013. 

48. Treissman D, Vidmar J, Guigard S, Kindzierski W, Schulz Jm Guigard E: Review 
of approaches for setting an objective for mixtures in ambient air using 
polycyclic aromatic hydrocarbons (PAHs). Edmonton, AB: Government of 
Alberta; 2004. 

49. Butler JP, Post GB, Lioy PJ, Waldman JM, Greenberg A: Assessment of 
carcinogenic risk from personal exposure to benzo(a)pyrene in the Total 
Human Environmental Exposure Study (THEES). Air Waste 1993, 
43:970-977. 



doi:1 0.1 1 86/1 476-069X-1 3-7 

Cite this article as: Irvine et al: Cancer risk to First Nations' people from 
exposure to polycyclic aromatic hydrocarbons near in-situ bitumen 
extraction in Cold Lake, Alberta. Environmental Health 2014 13:7. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at f~\ RiftMM i rpntral 

www.biomedcentral.com/submit momea central 



